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Abstract 

A persistent current, coupled with the spin state, of purely many- 
body origin is shown to exist in Nagaoka's ferromagnetic state in two 
dimensions (2D). This we regard as a manifestation of a gauge field, 
which comes from the surrounding spin configuration and acts on the 
hole motion, being coupled to the Aharonov-Bohm flux. This provides 
an example where the electron-electron interaction exerts a profound 
effect involving the spins in clean two-dimensional lattice systems in 

sharp contrast to continuum or spinless fermion systems. 
PACS numbers: 72.10.Bg, 05.30.Fk, 73.20.Dx 
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Persistent current is a highlight of the physics of electrons on mesoscopic scales. One 
major interest is how the electron-electron interaction affects the persistent current. The 
issue was originally posed to answer whether the interaction can offset the reduction of 
the current due to disorder. The many-body effect has been examined for both continuum 
systems and lattice systems with various geometries such as a one-dimensional ring or a 
torus i-i. 

In clean continuum systems, however, the electron-electron interaction, by itself, exerts 
null effect on the persistent current, which is exactly shown in one dimension (ID) and 
expected to persist in two dimensions (2D) and higher ones, unless the interaction creates 
low-energy excitations below the lowest excitation of the noninteracting system [^. 

The lattice system is a different story, even when clean, for two reasons. First, a lattice 
system can undergo Mott's metal-insulator transition at some band filling or at some critical 
interaction strength, so that the interaction may suppress the persistent current around the 
transition [p!0|,^. 

However, a truly interesting effect of interaction in lattice systems should appear when 
one considers the spin degrees of freedom. An effect of spin is known to appear as a change in 
the Aharonov-Bohm (AB) period (in units of the flux quantum, $o = hc/e) and amplitude, 
originally shown for the noninteracting model by averaging over the number of electrons 
||TT| and later found in interacting systems such as the one-dimensional Hubbard model 



with a fixed number of electrons in the dilute and/or strong-coupling limit. ||T2| Change in 



the orbital magnetism and current is also found in the finite-?/ Hubbard model with the 



spin-orbit interaction 13 



In the present paper we show that an interplay between charge and spin degrees of 
freedom can result in a qualitatively different phenomenon, i.e., a persistent current "driven 
by" a twist in the spin structure with the external flux in a strong-t/ Hubbard model in 2D. 

In strongly correlated systems, the ground state is often sensitive to the boundary con- 
dition. For instance, Nagaoka's ferromagnetic ground state which is an extreme man- 



ifestation of the correlation for large U in the immediate (one-hole doped) vicinity of the 



Mott insulator, evolves into a spin-singlet state when the antiperiodic boundary condition is 



adopted in one direction. We can then expect that the boundary- condition dependence 
will signify a large response to the AB flux and a persistent current via the Byers-Yang 
theorem. [|T^ This does indeed turn out to be the case. We shall trace this back to a 
gauge potential inherent in lattice systems, where the gauge potential has an ability to cope 
with the frustration caused by the twisted boundary condition, a situation distinct from 
continuum systems. 

We consider the two-dimensional Hubbard model on a torus to get rid of edges, and 
pierce an AB flux $ through the opening. This amounts to a frustration, or a change in the 
boundary condition along y for the wave function to 

0(. . . , Xi + NyViy, ...)= exp (i27r$/$o)0(- • • , x^, . . .) , 

where Xj is the position of the ith electron, and {No) the unit vector (the number of sites) 
in the a (= a;, y) direction with the lattice constant being unity. 

We have numerically obtained the change in the energy spectrum with $ at f/ = oo for 
15 (14) electrons on a 4 x 4-site system, i.e., one (two) -hole doped from half fllling) in Fig. 
|T]. We immediately recognize that there is a dramatic series of level crossings for low-lying 
states in such a way that the total spin of the ground state S oscillates between Nagaoka's 
full polarization (which occurs at $ = for one hole or at the half-flux for two holes) and 
vanishing (5* ^ 0) polarization. 

This oscillation gives rise to a persistent current / (Fig. 0) via the Byers-Yang theorem, 

m 

(9$ ' 

where |0) is assumed to be normalized. Thus, not only the magnitude but the overall proflle 
of the current deviate from the spinless electron system (the curve for the fully polarized 
state). 

The flux-induced oscillation in S found here is totally incompatible with the ordinary 
Stoner ferromagnetism. In the Stoner picture, the lowest charge excitation with a spin flip is 



the Stoner excitation, for which a Stoner gap opens over the whole Brillouin zone for a fully 
polarized state, so that the Stoner excitation cannot contribute to the persistent current. 
The collective spin wave also remains in an excited state when $ is varied. For partially 
polarized states, the Stoner excitation continuum becomes gapless and a finite $ can induce 
single spin flips, but this will only result in a change of S by unity. In Nagaoka's state, by 
contrast, S oscillates between the largest and smallest possible values, which is due to an 
anomaly in the gap and spin stiffness where the spin mass vanishes like the inverse system 

size MMM- 

In a doped Mott insulator, which we have at hand, a fundamental deviation from the 
spinless fermions arises from an interplay of the charge (hole) motion and the surrounding 
spin configuration, which is indeed responsible for anomalous excitations. To show this 
we can exploit the formulation of Ref. [0, where we write down an eigenfunction of the 
one-hole system in a "Bloch's form" as 

= Y: exp (ik ■ Xh) E /(^i, • • • , ^nJ 4,., ■ ■ ■ c^,.,^ |0). (2) 

Here Xh is the position of the hole, c]^^ creates an electron with spin aj at x^, and A'^e 
(= — 1) is the number of electrons in an A^-site system. The advantage of working with a 
spin wavefunction /(ai, . . . , ajyj around the hole is the following. If the order (xi, . . . , x^rj 
according to which the creation operators are multiplied, is specified by the coordinate 
relative to Xh, / no longer contains Xh, and the equation determining the spin configuration 
around the hole is simplified: A hole hopping from Xh to Xh + Uq is mapped onto a spin 
operation that is almost a uniform translation of the spins in the inverse direction, except 
for the spin at Xh + u^, which has to move to x^ — u^. If we write the whole operation as 
Ta, we end up with a spin Hamiltonian, 

d 

Kpin = -t ■Ua)Ta + CXp {-ik ■ U„) T^'^} . (3) 

a=l 

In this formalism, $ simply shifts the k points uniformly as k ^ k + (0, 27r$/(A^y$o))- 
We can then identify the origin of the many-body effect on the current as the gauge 
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potential in the following sense. When we decompose the motion of electrons into the center- 
of-mass and internal ones, $ only couples to the former, which is completely decoupled from 
the latter for continuum systems. Thus, a persistent current cannot deviate from the free- 
fermion behavior 0. Even if the interaction creates low-energy excitations leading to level 
crossings for finite the curvature of the spectral flow does not deviate from that 

for the free fermions. 

In sharp contrast, the center-of-mass motion does disturb the relative coordinates on 
a lattice. Namely, a shift in Xh, which coincides with the center-of-mass here, affects the 
relative coordinates between electrons through the electron correlation (i.e., excluded double 
occupancies). In the present formalism, the coupling of the two coordinates appears as a 
discrepancy of Ta from the uniform translation. |[T^ The spin wave function is thus affected 
by the fiux-coupled charge motion, and changes its structure with $ in two- (or higher) 
dimensional lattice systems. 

Here we notice that the spin configuration does not flow into an antiferromagnetic (AF) 
state with $, since this would impede the motion of the hole. In fact, the state that takes over 
Nagaoka's state at $o/2 is a twisted spin structure, or a spin-density wave whose wavelength 
is as large as the sample size. The spiral structure is confirmed here from the spin-spin 
correlation for a (v^26 x ^/26)-site system displayed as an inset of Fig. |l](a), where the spin 
polarization is seen to rotate gradually along the y direction in which the AB flux exerts a 
twist. 

Curiously, it is exactly this spiral state that is the first excited state from Nagaoka's at 
$ = 0. Among the lowest-spin {S = 1/2) states, this (transfer-stabilized) spiral-spin state 
has a lower energy than the (exchange-stabilized) AF state until a level crossing occurs at a 
critical value of U [(= 167t, i.e., J = At'^/U = 0.024t for a (^26 x y26)-site system)] as U 
is decreased. 

The same occurs for a two-hole system, where the ground state starts off with the spiral 
spin state at $ = (for small enough J). The state is characterized by a spin-spin 
correlation S'(Q) having four peaks at the wave number Q = (±27r/A^^., ±27r/A^j^) [inset of 
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Fig. 1^ (b)]. The 5* = state changes into a two-peaked one continuously with $. 

Such twisted spin structures have in fact been analyzed to discuss the stability of Na- 
gaoka's state. Dougot and Wen have introduced a twisted spin state, where they assume 
that a hole hops on a frozen spin background. There, they have shown that a properly cho- 
sen spin texture can remove some of the frustration caused by the external flux. Although 
the Dougot-Wen function is an approximate trial one, which is not even an eigenstate of 
the total spin, the function, when cast into our Bloch form, does serve a heuristic purpose 
of identifying two kinds of Bloch momenta. We can first call k, which is coupled to Xh in 
Eq. (H), the charge Bloch wave number. The Dougot-Wen spiral configuration is frozen in 
the rest frame, but is now specified by the position of its loop fq = (xc^/o) (where the spin 
wave function has a phase slip of vr) in our frame moving with the hole. We can take a 
linear combination of spin configurations having different tq with a spin Bloch wave number 
q. Let 6 = 2TT/Ny he the wave number of the spin twist in the Dougot-Wen function. The 
variational function then reads 

cos[f (y + yo)] \ 



|^)=Eexp(ik-x,)^Ee^''-^°n' 



(4) 

sin[|(?/ + ?/o)] / (^'f) 



where \'^)(x,y) is the spinor at site {x,y) and we have again assumed that the order in the 
product Y[{x,y) is according to the rule described below Eq. (0). If we choose q = (0,9/2) 
for the spin Bloch wave number, the variational energy becomes 

e /, 2tt^ e\ 

e = -2t cos kx - 2t cos - cos \ ky + I , (5) 

for U = oo, and is minimized at the half-fiux, where the difference between the AB phase 
271^/Ny^Q and the spin phase 6/2 vanishes. 

Thus we are led to a picture in which the extra phase 6/2 arising from the gauge potential 
(an internal frustration) from the spin configuration acts to cancel the enhancement in the 
kinetic energy caused by the external (fiux-induced) frustration. There, the low-spin states 
becomes advantageous for the maximal (half-fiux) frustration, since they have the maximum 
degrees of freedom in the spin space that enable them to absorb the frustration. 
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We expect this should generally apply to strongly correlated systems with restricted 
motion of charges (i.e., excluded double occupancies), so that the variation of the energy, 
and thus the persistent current, for spinfuU electrons should be smaller than those for (the 
same number of) spinless electrons. 

To be more precise, the spin configuration is disturbed by the hole motion as stressed 



2^ , and a factor neglected by the Dougot-Wen wave function is the quantum fluctuation. 
This is numerically seen in a reduction in the absolute value of the nearest-neighbor spin-spin 
correlation from that for the spiral configuration of classical S = 1/2 spin in the inset of Fig. 
|l](a). However, the inseparability of the center-of-mass and internal coordinates and their 
interference does in general remain when the configurational subspace classified by the total 
orbital angular momentum and the spin are 'irreducible' (i.e., when any two states within 
the subspace are connected by a series of nonzero matrix elements of the transfer). Thus, in 
addition to the well-known enhancement in the effective mass of a hole, this gives rise to the 
spin-modulated persistent current in clean two-dimensional lattices, unlike two-dimensional 
continuum systems. 

In other words, a current deviating from the free-electron behavior should be generic to 
any lattice systems, such as the t-J model, that have a constrained electron motion due to 
strong interactions. For example, in the AF phase in the finite-t/ Hubbard model or t-J 



model in 2D, £'»($) deviates from the free-electron result as shown by Poilblanc [^. A 
difference there from the present situation is that the level crossing is shown to occur among 
the exchange-stabilized Stot = states (i.e., within the AF phase). There, we also observe a 
reduction in the deviation of the ground-state energy, which may again be explained from 
the flux-induced enhancement of the kinetic energy as relaxed by the gauge field arising from 
the spin wave function. 

Obviously, another important question is the effect of disorder or the edge effect for 
open geometries to make the model approach to actual mesoscopic systems. As for the 
ferromagnetism, whether there is a finite ferromagnetic region around Nagaoka's limit is 
still an open question. A finite region for the partial polarization is proposed from a t- 



J study but a possibility of an abrupt transition to an AF state is also suggested 
pO| . However, we can provide a way to relax the singular situation of Nagaoka's by going 
over to two-band Hubbard models, where the double- exchange mechanism induces a more 



stable ferromagnetism that we can show has close similarities with Nagaoka's . Another 
prominent example of the stabilized ferromagnetism is the Hubbard model on fiat-band 
systems |26| . There the fiat band, which is fully polarized when half-filled, becomes metallic 
when doped, and it will be interesting to see how the persistent current behaves in such an 
itinerant magnet. 
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FIGURES 

FIG. 1. (a) AB flux ($) dependence of the lowest state for each of various total spins in 
a one-hole doped (4 x 4)-site system with U = oo, in which the ground state evolves from a 
ferromagnetic state at $ = down to a spiral-spin state at $ = $o/2. (b) The same for a two-hole 
doped (4 X 4)-site system. The insets depict the spin-spin correlation function {S{0)S{'R)) for the 
spiral-spin ground state with the half-flux in a one-hole doped (\/26 x \/26)-site system [the arrow 
in (a)] and the spin structure factor S'(Q) for the ground state without flux in a two-hole doped 
(\/20 X \/20)-site system [the arrow in (b)] at C/ = oo (J = 0). 

FIG. 2. The persistent current derived from the previous figure via the Byers-Yang theorem 
for the (4 x 4)-site Hubbard model for one hole (a) or two holes (b) with U = oo. 



11 



